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We investigate an inequality constraining the energy and potential enstrophy flux in the 
two-layer quasi-geostrophic model. This flux inequality is unconditionally satisfied for 
the case of two-dimensional Navier-Stokes turbulence. However, it is not obvious that it 
remains valid under the multi-layer quasi-geostrophic model. The physical significance of 
this inequality is that it decides whether any given model can reproduce the Nastrom- 
Gage spectrum of the atmosphere, at least in terms of the total energy spectrum. We 
derive the general form of the energy and potential enstrophy dissipation rate spectra for 
a generalized multi-layer model. We then specialize these results for the case of the two- 
layer quasi-geostrophic model under dissipation configurations in which the dissipation 
terms for each layer are dependent only on the streamfunction or potential vorticity of 
that layer. We derive sufficient conditions for satisfying the flux inequality and discuss 
the possibility of violating it under different conditions. 



1. Introduction 

It is now well-known that in two-dimensional Navier-Stokes turbulence, most of the 
energy tends to go towards larges scales and most of the enstrophy tends to go towards 
small scales. This was initially proposed by Fj0rt0ft (19531 via his triad interactions 
argument. Later, Kraichnan (1967), Leith (1968), and Batchelor (1969) introduced the 
theory that the energy forms an upscale inverse energy cascade with energy spectrum 
scaling as fc _5//3 , and that the enstrophy forms a downscale enstrophy cascade with k~ 3 
scaling. Kraichnan (1967) argued, differently from Fj0rt0ft (1953), that the direction of 
the two cascades can be justified via a thermodynamic argument in which we introduce, 
without proof, the assumption that the energy and enstrophy fluxes should tend to 
revert the energy spectrum from a cascade configuration to the absolute equilibrium 
configuration. The existence of forcing and dissipation arrests this tendency, thus keeping 
the system locked in a steady-state forced-dissipativc configuration away from absolute 
equilibrium. 



Merilees & Warn ( 1975 ) identified a serious error with the original Fj0rt0ft argument: 



Fj0rt0ft (1953) claimed that in every triad interaction group, more energy is transferred 



upscale than downscale. However, a more rigorous analysis shows that there exist triad 
interaction groups in which more energy is sent downscale than upscale, and it is not 
obvious, without additional considerations, as to which group is dominant. This was 
explained in detail also by Gkioulekas & Tung (2006). Aside from this matter, the funda- 
mental problem that underlies every other proof that utilizes only the twin conservation 
laws of enstrophy and energy, is that an additional assumption needs to be introduced to 
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overcome the symmetry of the Euler equations under time reversal. Typical assumptions, 
such as the tendency of the energy spectrum to revert to absolute equilibrium, or the 
tendency of an energy peak to spread, typify ad hoc constraints imposed implicitly on 
the initial conditions that are needed to break the time reversal symmetry. All of these 



proofs have been reviewed in detail by Gkioulekas & Tung (20076). More importantly, 



Gkioulekas & Tung (20076 ) have counterproposed a very simple and mathematically rig- 



orous proof that avoids the need for any ad hoc assumptions by considering the combined 
effect of the Navier-Stokes nonlinearity and the dissipation terms. The only assumption 
used by this proof is that the forcing spectrum is restricted to a finite interval [fcijfe] 



of wavenumbers. Later, Farazmand (2010) and Farazmand, Kevlahan & Protas (20111 



showed that even that assumption can be relaxed to some extent, although not entirely 
eliminated. 

The essence of the argument is to show that for every wavenumber k not in the forc- 
ing range, the energy flux TlE(k) and the enstrophy flux Ha(k) satisfy the inequality 
fc 2 n^;(fc) — IlGf(fc) < 0. Here, IT£(fc) represents the amount of energy per unit vol- 
ume transferred from the wavenumbers in the (0,fc) interval to the wavenumbers in the 
(k, +oo) interval, and Hc(k) is defined similarly for the enstrophy. From this inequality 
we then derive the following integral constraints for TlE(k) and Hc{k): 



[ qn E (q)dq<0,\fke{k 2 ,+oo), 
Jo 

r+oo 

/ q- 3 U G (q)dq>0, Vk £ (0, fc X ). 
Jk 



(1.1) 

(1.2) 



As explained by Gkioulekas fc Tung| (2007&), these constraints imply a predominantly 
upscale transfer of energy and a predominantly downscale transfer of enstrophy. The 
original flux inequality k 2 HE(k) — Hc(k) < itself can also be directly interpreted as a 
tight constraint on the downscale energy flux. 

The question that concerns us in the present paper is whether this flux inequality con- 
tinues to hold for multi-layer quasi-geostrophic models. Quasi-geostrophic turbulence, as 



shown recently by Lindborg (2007), is a reasonable approximation for the cascade dynam- 



ics of atmospheric turbulence at length-scales larger than 100km, in order of magnitude. 
Underlying the model arc the assumptions of small thickness and fast rotation. Towards 
large scales, the atmospheric thickness becomes relatively small and the Coriolis effect 
of rotation becomes increasingly noticeable. In multi-layer quasi-geostrophic models, the 
vertical dimension is discretized into a finite number of layers. The simplest layer model 
is one in which we only have two layers. This two-layer quasi-geostrophic model was simu- 



lated by Tung & Orlando ( 2003 ), who claimed that it can reproduce the energy spectrum 



of the atmosphere, consistently with the measurement of Gage & Nastrom (19861. Tung 



& Orlando (2003) intended to corroborate their theory that the Nastrom-Gage spectrum 



consists of a downscale potential enstrophy cascade coexisting with a downscale energy 
cascade over the same range of scales. Gkioulekas & Tung (2005a b ) clarified the theoret- 
ical mechanism that allows cascades to coexist and predicted that each cascade provides 
a contributing term to the energy spectrum and that both contributions are linearly su- 
perposed with one another. This linear superposition hypothesis and the Tung-Orlando 
theory of coexisting cascades have been both corroborated by the recent measurements 



and analysis by Terasaki, Tanaka & Zagar (2011 ). A recent primitive equations numerical 



simulation by Vallgren, Deusebio & Lindborg (2011) also seems to support the notion 



of a downscale energy cascade coexisting with a downscale potential enstrophy cascade, 



and a more detailed discussion of these matters has been given by Gkioulekas (2012). 
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More specifically, the Nastrom Gage spectrum of the atmosphere consists of a k~ 3 
slope in the energy spectrum that transitions to k~ 5 ^ 3 at a length scale of approximately 
1000km to 700km, near the Rossby deformation wavenumber k^. According to the theory 
of |Tung fc Orlando| ( [2003| and |Gkioulekas fc Tung| ( |2005a|p| ), the Nastrom-Gage spectrum 
arises from coexisting downscale cascades of potential enstrophy and energy, over the 
same range of scales, each cascade gives a contribution to the energy spectrum, and both 
contributions are combined linearly. As a result, the k~ 3 contribution of the downscale 
potential enstrophy cascade is dominant for small wave numbers fc, only to be overtaken 
by the fc -5 / 3 contribution of the downscale energy cascade near a transition wavenumber 
fc t . For the fc~ 5 / 3 contribution of the downscale energy cascade to become dominant 
at wavenumbers fc ~ fc t beyond a transition wavenumber fc t , it is necessary for the 
flux inequality to break down near kt and reverse direction for wavenumbers k > k t . 
Consequently, the direction of the flux inequality decides whether various models can 
reproduce the Nastrom-Gage spectrum, at least in terms of the total energy. 

For the case of two-dimensional turbulence, as was first noted by Gkioulckas & Tung 



( 2005a|6 1, the unconditional validity of the flux inequality constitutes a matheniati 
cally rigorous proof that two-dimensional turbulence cannot reproduce the Nastrom-Gage 
spectrum under any circumstances, as long as the dissipation terms are given in terms of 
integer or fractional powers of Laplacian operators. For the case of the real atmosphere, it 
is fairly obvious that the flux inequality is violated, provided that we agree with the Tung- 
Orlando theory ( Tung fc Orlando|2003 ) that the Nastrom-Gage spectrum consists of two 
coexisting downscale cascades of potential enstrophy and energy. If the flux inequality 
can be also violated by quasi-geostrophic models, then these models can reproduce the 
Nastrom-Gage energy spectrum, at least in terms of the total energy spectrum E(k), if 
not in terms of a correct distribution of energy between kinetic energy and potential en- 
ergy. If, on the other hand, the flux inequality continues to hold, then quasi-geostrophic 
models cannot reproduce the Nastrom-Gage spectrum either. The question that we ulti- 
mately want to investigate is whether the quasi-geostrophic models have enough physics 
to enable a break-down of the flux inequality, and what configuration of the dissipation 
terms is needed to accomplish that. 



Gkioulekas & Tung (2007a) first noticed that the only mechanism that can cause 
a flux inequality violation is an asymmetry in the dissipation operators for different 
layers. Unfortunately, without more information about the phenomenology of the two- 
layer quasi-geostrophic model, and multi-layer quasi-geostrophic models in general, it is 
very difficult to derive any sufficient conditions for the needed flux inequality breakdown. 
On the other hand, it is possible to derive rigorous sufficient conditions for satisfying the 
flux inequality when the asymmetry is restricted by an upper bound, without assuming 
any knowledge whatsoever of the phenomenology of any of the spectra involved. The goal 
of the present paper is to address this non-controversial aspect of the larger problem, by 
deriving a series of such results, for the case of weak dissipation asymmetry, and in doing 
so, to also set up the necessary mathematical infrastructure for a further exploration of 
the flux inequality. Unlike the case of two-dimensional turbulence, where the derivation 
of the flux inequality is easy, it should be clear from the present manuscript that the 
flux inequality in the quasi-geostrophic model is both nontrivial and sensitive to the 
configuration of the dissipation terms across layers. 

This paper is organized as follows. In section 2, we introduce the generalized multi- 
layer model and its conservation laws of energy and layer-by-layer potential enstrophy. 
In section 3, we introduce the budget equations for the energy and potential enstrophy 
spectra, derive the required conditions for the energy spectrum definition to remain 
positive definite, and derive general expressions for the energy and potential enstrophy 
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dissipation rate spectra. Simplifications of these general expressions are given for the 
special cases of streamfunction dissipation and potential vorticity dissipation. From these 
results, in section 4 we derive sufficient conditions for satisfying the flux inequality for 
the special case of the two-layer quasi-geostrophic model under the two special cases 
of streamfunction dissipation and potential vorticity dissipation. Because the details of 
the calculation are tedious, the main results are summarized in section 5. Conclusions 
and general discussion are given in section 6. Technical matters are taken up in the 
appendices. 



2. The generalized multilayer model and conservation laws 



Following Gkioulekas (|2012|), we write the governing equations for the generalized 

+ J(lp a ,q a ) =d a + fa, 



multi-layer model in matrix form: 

dq a 
dt 



(2.1) 
(2.2) 



Here ip a represents the streamfunction at the a-layer, q a represents the potential vorticity 
at the a-layer, is a linear operator encapsulating the dissipation terms, and f a is 
the forcing term acting on the a-layer. The index a takes the values a = l,2,...,n 
representing the layer number, for a model involving n layers. Sums over indices, such as 
in the sum over the index j3 in the dissipation terms above, are assumed to run over all 
layers 1,2, ... ,n, unless we indicate otherwise. It is also assumed that the streamfunction 
ij) a and the potential vorticity q a are related via a linear operator _£? a| g according to: 



9a(x,t) 







^a/3^a (x,i). 



(2.3) 



The above equations encompass both the two-layer quasi-geostrophic model and the mul- 
tilayer quasi-geostrophic model, on the assumption that we neglect the /3-effect, arising 
from the latitudinal dependence of the Coriolis pseudoforce. This is a reasonable assump- 
tion for Earth, especially if we restrict our interest to a thin strip of the Earth's surface, 
oriented parallel to the equator. Baroclinicity instability is accounted for by the forcing 
term f a , and implicit in the entire argument is the assumption that it forces the system 



at large scales only. This assumption, originally proposed by Salmon (1978 1980), is the 



only physical assumption implicit in the theoretical framework of the flux inequality, and 



it has been corroborated by Welch fc Tung (1998) and Tung fc Orlando (20031 



For the sake of simplifying our analysis, we assume that all fields are defined in an infi- 
nite two-dimensional domain. Then we can write the Fourier expansions for the stream- 
function ip and the potential vorticity q as follows: 



V> Q (k, t) exp(ik ■ x) dk, 



q a (x,t) = / <? a (k, t) exp(ik • x) dk. 

Jm. 2 



(2.4) 
(2.5) 



We assume that the operator Jzf Q/ 3 is diagonal in Fourier space. This means that the 
relation between the streamfunction and the potential vorticity, in Fourier space, reads: 



q a (k,t) 







|k||hMM). 



(2.6) 
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We also assume that J£ a p is symmetric with Jzf a ^ = J£p a . This implies that L a p(k) = 
L/3 a (k) for all wavenumbers k. For an n-layer quasi-geostrophic model, the L a p(k) matrix 
becomes a tridiagonal matrix. For the special case of the two-layer quasi-geostrophic 
model, the matrix L a p(k) is given by 



L a p(k) = 



-k' 2 



kl/2 



+4/2 



i/2 



(2.7) 



with fcfl the Rossby wavenumber. For quasi-geostrophic models, the matrix L a p{k) is 
non-singular for all wavenumbers k > 0, due to being diagonally dominant, and we 
assume that to be the case in our abstract formulation given above. Consequently, there 
is an inverse matrix L~^(k) which defines the inverse operator J/f~p . To accommodate 
a possible singularity at k — we assume that at wavenumber k = 0, in Fourier space, 
the corresponding field component is for all fields. This is equivalent to subtracting the 
mean field and considering only the field fluctuation around the mean. 

Likewise, we assume that the dissipation operation %} a p is also diagonal in Fourier 
space and that the Fourier expansion of the dissipation term 3> a /3ip/3 reads: 



-D a/ 3(||k||)Vi / 3(k,t)exp(ik • x) dk. 



(2.8) 



In general, the dissipation operator is a matrix operator that is allowed to entangle 



multiple layers together. Realistically, this occurs when one follows Salmon ( 1978 1980 ) 



in the definition of the Ekman term, which should appear only on the bottom layer 



but entangles both layers together. Previously, Tung & Orlando (2003), Gkioulekas & 
Tung (2007a), and Gkioulekas ( |2012 ), restricted their attention to diagonal dissipation 
only where D a p{k) = for a ^ /3. In this case, the spectrum D a p(k) of the dissipation 
operator reads D a p{k) = 5 a pDp{k) with 8 a fj given by 



1, ifa = P 
0, if a ^ /3 



(2.9) 



We will now show that the generalized layer model, in the absence of dissipation, 
conserves the total energy E and the total potential enstrophy G under very general 
conditions on the operator Jzf a /3, For any arbitrary scalar field f(x,y) we write the cor- 
responding volume integral using the following notation: 



((f)) 



f(x,y) dxdy. 



(2.10) 



We define the total energy E over all layers, and the layer-by-layer total potential en- 
strophy G a for layer a, as E = — ^a^aQa)) and G a — ((c?^)). The purpose of the 
minus sign in our definition of the total energy E is to maintain consistency with the 



notation and sign conventions used by Gkioulekas (2012). Specifically, we will show that 



the potential enstrophy is conserved on a layer-by-layer basis unconditionally regardless 
of the details of the operator Jz? Q ^. Conservation of the total energy E, over all layers, 
on the other hand, requires that the operator ££ a fi be symmetric and self- adjoint. By 
symmetric we mean that the operator satisfies J? a p = J£p a . To define the self-adjoint 
property, consider two arbitrary two-dimensional scalar fields f(x, y) and g(x, y). We re- 
quire that every component of the operator J£ a p must satisfy ((f(^f a pg))) = (((Jf a pf)g)) 
for any two fields f(x,y) and g{x,y). This self-adjoint property, so defined, follows as an 
immediate consequence of our previous assumption that the operator ££ a p is diagonal in 
Fourier space. In the proof given below, however, there is no need to use the stronger 
assumption of diagonality. 
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The proof is based on the following properties of the nonlinear Jacobian term. If a(x, y) 
and b(x, y) are two-dimensional scalar fields that satisfy a homogeneous (Dirichlet or Neu- 
mann) boundary condition, then we can show that ((J(o,6))) = 0, using integration by 
parts. Then, we note that, as an immediate consequence of the product rule of differen- 
tiation, given three two-dimensional scalar fields a(x,y),b(x 7 y), and c(x,y) we have 

((J(ab, c)» = ((aJ(b, c))) + ((bJ(a, c))> = 0, (2.11) 

from which we obtain the identity 

«aJ(&,c)» = ({bJ(c,a}}) = ({cJ(a,b)}). (2.12) 

Now, let us go ahead and drop the dissipation and forcing terms and write the time- 
derivative of the potential vorticity q a as q a = —J(ijj a ,q a )- Then, the time derivative of 
the streamfunction ip a reads: 

i>* = E = - E #$Wp> <&» )• ( 2 - 13 ) 

Differentiating the total potential enstrophy G a for the a layer with respect to time 
and employing the identity given by Eq. (2.12) immediately gives: 

G« = 2((q a q a )) = -2{{q a J{ij} a ,q a ))) = -2{{i/> a J(q a ,q a ))) = 0. (2.14) 

Here, we note that from the definition of the Jacobian J(q a ,q a ) — 0. This establishes 
the layer-by-layer conservation law of potential enstrophy, unconditionally, as claimed. 
To show the energy conservation law, we differentiate the total energy E with respect to 
time and obtain: 

E = -(d/dt) X><*fc» = - £)«lMa» - £)«lMa» (2.15) 
a a a 

= J2((l a ^J^^qp)))+J2^ J (^q a ))) (2.16) 

af) a 

= 52{{JU>p,qf})-ZZpq a )) + J(i>a^a))) (2.17) 

af) a 

= J2{(J(iPf,,qp)J?^q a )) = J2((J(^QP)^)) (2-18) 

a/3 

= Y / ((J(Mp)Q0)) = O- (2.19) 



Note that the self-adjoint property is applied at step Eq. (2.17), and the symmetric 
property is applied at Eq. (2.18). This concludes the proof. 



3. Spectra and budget equations 



Following Frisch (19951 and Gkioulekas (20121, we define spectra for the energy and 



potential enstrophy using the bracket notation, which is defined as follows. Consider, 
in general, two arbitrary two-dimensional scalar fields a(x) and b(x). Let a <fc (x) and 
b <k (x) be the fields obtained from a(x) and 6(x) by setting to zero, in Fourier space, 
the components corresponding to wavenumbers whose norm is greater than k. Formally, 
a <fe (x) is defined as 



a <fc (x) 



dx / dk 



g(fc-||k Q ||) 
4tt 2 



exp(ik • (x - x ))a(x ) 



(3-1) 
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with H(x) the Heaviside function, defined as the integral of a delta function: 

fX f 1, if if x G (0, +00) 

H(x) = / <J(t) dr = ) 1/2, if if x = . (3.2) 

Jo { 0, if if a; G (-00, 0) 

Obviously, 6 <fc (x) is defined similarly. We now use the two filtered fields a <fe (x) and 
6 <fc (x) to define the bracket (a,b) k as: 



(a,b) 



k dk 
1 



dx (a <k (x)b <k (x.)) 

dCl(A) ([a* (kAe)b(kAe) + a{kAe)b*{kAe)} 



Aeso(2) 



(3.3) 
(3.4) 



Here, d(k) and 6(k) are the Fourier transforms of a(x) and 6(x), SO(2) is the set of all 
non-reflecting rotation matrices in two dimensions, dCl(A) is the measure of a spherical 
integral, e is a two-dimensional unit vector, and (•) represents taking an ensemble average. 



The star superscript represents taking the complex conjugate. Note that Eq. (3.3) is the 



definition of the bracket, and Eq. (3.4) follows from Eq. (3.3) as a consequence. As noted 



by Gkioulekas (2012), the bracket satisfies the following properties: 

(a,b) k = (b,a) k , 

(a,b + c) k = (a,b) k + {a,c) k , 

(a + b,c) k = (a,c) k + (b,c) k . 



(3.5) 
(3.6) 
(3.7) 



Moreover, every (a/3)-component of the operator Jz? Q( 3 is self-adjoint with respect to the 
bracket, which gives 



xpa,b) k = (a,J!? a pb) k = L a/3 (k) (a,b) k , 



(3.8) 



and the same property is also satisfied by every component of the inverse operator J2? a } : 



Kp^) k = ( a >^P b ) k = L ^ k ) M) fc ■ 



(3.9) 



Using the bracket, we define the energy spectrum E(k) = — {ip a , <Za)fcj and we also 
define the layer-by-layer potential enstrophy spectrum G a (k) — (q a ,<la} k and the total 
potential enstrophy spectrum G(k) — ^2 a G a (k). Unlike the case of two-dimensional 
Navier-Stokes, where the enstrophy and energy spectra G(k) and E(k) are related via a 
simple equation, G(k) = k 2 E(k), in the generalized layer model, the potential enstrophy 
spectrum and the energy spectrum are related indirectly, as shown below: 

Define the streamfunction spectrum C a p(k) = (ip a ,il>p} k . Then, via the properties of 
the bracket above, the energy spectrum E(k) reads 



a a \ / 

= -);L a p(k)C a p(k), 

a/3 

and the potential enstrophy spectrum G a (k) reads 



a a \ P 7 I k 



Y,L a p{k){^ a ,^p) k (3.10) 

up 



(3.11) 



(3.12) 
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= L ot/3(k) ( i>/3, J~] ^7^7 ) = ^2L al3 (k)L aj (k) (ipp,ip-y) k (3.13) 

a/3 \ 7 / fc aPf 

= J2 L Mk)L ai (k)C^(k). (3.14) 

a/37 

Thus, they are related only indirectly via the streamfunction spectrum C a p(k). This was 



noted previously by Gkioulekas (2012). 



We note that for a ^ (3, C a p{k) is not positive-definite and may take positive or 
negative values. For the case a = (3 we define U a (k) = (ij} a , ip a )k, which is positive 
definite (i.e., U a (k) > 0). Then we note that 2\C a p(k)\ < U a (k) + Up(k). We can use this 
inequality to show that if the matrix L a p(k) satisfies the diagonal dominance condition 

L a p{k) > 0, for /3, (3.15) 

2L aa (k) + (L aP {k) + L Pa {k)) < 0, (3.16) 

a/3 
a^f) 

then the energy spectrum E{k) is positive definite. We give the proof in Appendix |A| 
Both the two-layer quasi-geostrophic model and the multi-layer quasi-geostrophic model 
satisfy this diagonal dominance condition. As for the layer-by-layer potential enstrophy 
spectra G a (k), it is immediately obvious that they are unconditionally positive definite, 
regardless of the form of the matrix L a p{k), since by definition G a {k) = (q a , q a ) k . 

For the argument of the present paper we need the dissipation rates D E {k) and Dc(k) 
for the energy and total potential enstrophy expressed in terms of the streamfunction 
spectrum C a f}(k). In appendix [B] we show that the energy dissipation rate spectrum 
DE(k) and the layer-by-layer potential enstrophy dissipation rate spectra Dc a (k) are 
given by 

D E (k) = 2j2 D ap(k)C a0 (k), (3.17) 

D Ga (k) = -2 Y, L a p{k)D an {k)C M {k). (3.18) 

Note that in order for the dissipation terms to be truly dissipative, the dissipation spec- 
tra Dsik) and Dc(k) need to be both positive definite. From the general form of the 
above equations this is not readily obvious. However, for simpler configurations of the 
dissipation operators, the above expressions for D E (k) and Da(k) simplify considerably, 
thereby making it possible to establish positive definitcness. 

We restrict our attention to cases where the dissipation operators at every layer involve 
only the fields of the corresponding layer, with no explicit interlayer terms. This can be 
arranged in terms of a linear operator 3> a applied to either the streamfunction i\) a or 
the potential vorticity q a . If D a {k) is the spectrum of the positive-definite operator 
then for the case of a dissipation term d a = £l a ip a , we have D a p(k) = S a pDp(k). We 
designate this case as streamfunction-dissipation. The D E {k) and Dc a {k) simplify as: 

D E {k) = 2j2D a p{.k)C a p(k) = 2Y,8 afi D p (k)C aP (k) (3.19) 

a/3 a0 

= 2Y,D a {k)C aa {k) = 2Y / D a (k)U a (k) 1 (3.20) 

a a 

D Ga (A!) = -2^£^(fc)I? Q7 (A)q^(*) = -25^L a/J (*)J crr £> 7 (A)C , /l7 (A). (3.21) 

M Pi 



the two-layer quasi- geostrophic model 
-2^L a0 (k)D (k)C a ^k). 



9 

(3.22) 



Note that for D a (k) > 0, it follows that D E {k) > 0, but it is not obvious that the same 
result extends to D Ga (k). However, if we further assume that the same operator is used 
for all layers, i.e. D a (k) = D(k), then we have the more specialized case of symmetric 
streamfunction- dissipation, and the dissipation rate spectra D E (k) and Dc{k) can be 
simplified further to give: 



D E (k) = 2J2 D a (k)U a (k) = 2D(k) J2 U a (k) = 2D(k)U(k), 

a a 

D G (k) = Y / D G Jk) = -2Y,L a p(k)Dp(k)C a p(k) 

a a/3 



2D(k) 



-^2L aP (k)C a p{k) 

a/3 



2D(k)E(k). 



(3.23) 
(3.24) 

(3.25) 



Now, D(k) > implies both D E {k) > and D G {k) > 0. 

Another possible arrangement is potential vorticity dissipation in which the positive- 
definite operator & a is applied to the potential vorticity q a , thereby yielding a dissipation 
term of the form 



d a = -S> a q a = - ^a^a/3^/3 = ^ ■ 



(3.26) 



It is easy to see that in this case, D a p(k) = —D a (k)L a p(k), and consequently, the 
dissipation rate spectra D E (k) and Dc(k) simplify as follows: 



D E {k) = 2j2D af 3(k)Caf3(k) = -2j2D a (k)L af} (k)C a p(k), 



(3.27) 



a/3 



D G Jk) = -2Y,L a p(k)D ai (k)C^(k) = -2j2L a0 (k)[-D a (k)L a7 (k)]C^(k) (3.28) 

f3j M 



2D a (k) 



Y L a /3{k)L ai (k)C 01 {k) 



2D a (k)G a (k). 



(3.29) 



Now, for D a (k) > 0, it follows that D Ga (k) > 0, whereas it is now the sign of D E (k) that 
is uncertain. Similarly to the previous case, if we assume symmetric potential vorticity 
dissipation, where the same dissipation operator is used on all layers, we then set D a (k) = 
D(k), and our expression for the dissipation rate spectra D E (k) and D G (k) simplify to 



D E (k) = -2j2D a (k)L afj (k)C aP (k) 

a/3 



= 2D(k) 



^ L a p(k)C a p(k) 



a/3 



= 2D{k)E{k), 



D G (k) = 2j2D a (k)G a (k) = 2D(k) £ G a (k) = 2D(k)G(k). 



(3.30) 

(3.31) 
(3.32) 



Now, for D{k) > 0, it follows again that both D E (k) > and D G (k) > 0. Note that the 
negative sign in —S$ a q a was necessary to ensure that the dissipation term d a = —S> a qa 
remains dissipative when the operator 3> a is positive-definite. 
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4. Flux Inequality for the two-layer model 

We now turn to the main issue of identifying sufficient conditions for satisfying the 
flux inequality k 2 HE(k) — Ha(k) < for quasi-geostrophic models. Presently, we re- 
strict our interest to the two-layer quasi-geostrophic model under symmetric or asym- 
metric streamfunction-dissipation or potential vorticity-dissipation. The multilayer case 
and other configurations for the dissipation terms will be considered in future publica- 
tions. We begin with a brief review of the two-dimensional Navier-Stokes case and the 
formal setup of the argument in section |4.1| Our calculations for the case of streamfunc- 
tion dissipation are given in section [4~2} and for the case of potential enstrophy dissipation 
in section |4.3| Since the details of the calculations are somewhat technical, a summary 
discussion of the main results is given in section [5] 



4.1. Preliminaries 

Let us recall that the energy flux spectrum TlE(k) is defined as the amount of energy 
transferred from the (0, k) interval to the (fc, +oo) interval per unit time and per unit 
volume. Likewise, the potential enstrophy flux spectrum Tl G (k) is the amount of potential 
enstrophy transferred from the (0, k) interval to the (k, +oo) interval, again per unit time 
and volume. Assuming a forced-dissipative configuration at steady state and that the 
wavenumber k is not in the forcing range, the energy and potential enstrophy transferred 
into the (fc, +oo) interval eventually are dissipated somewhere in that interval. It follows 
that we may write the flux spectra Hb(£;) and Hoik) as integrals of the energy and 
potential enstrophy dissipation rate spectra -D £;(£;) and D G (k): 

r+oo 

U E (k)= / D E (q)dq, (4.1) 
Jk 



k 



lI G (fe) = / D G (q)dq, (4.2) 

which implies that 



/•-too p+oo 

k 2 U E (k)^U G (k) = [k 2 D E (q)-D G (q)]dq = A(k, q)dq. (4.3) 

Jk Jk 

We see that a sufficient condition for establishing the flux inequality is to show that 
A(k,q) < for all wavenumbers k < q. It is also easy to see that A(fc, q) > for all 
wavenumbers k t < k < q is sufficient for establishing the violation of the flux inequality 
for all wavenumbers k > k t . 

For the case of two-dimensional Navier-Stokes turbulence, the dissipation rate spectra 
DE(k) and D G (k) are related via D G (k) = k 2 DE(k). This immediately gives A(k,q) = 
k 2 D E (q) — D G (q) = (k 2 — q 2 )D E (q) < for all wavenumbers k < q (since D E (k) > 0), 
which in turn gives the flux inequality k 2 HE(k) — H G (k) < 0. The physical interpretation 
of this inequality is that when we stretch the separation of scales in the downscale range, 
the energy dissipation rate at small-scales vanishes rapidly. As a result, most of the 
injected energy cannot cascade downscale although, as noted by |Gkioulekas fc Tung] 
(200 5a|6 1, a small amount of energy is able to do so. As we have seen in the previous 



section, for the case of quasi-geostrophic models, the energy and potential enstrophy 
dissipation rate spectra no longer have a direct and simple relation with each other, so 
the validity of the flux inequality needs to be carefully re-examined. 

We set up our analysis of the two-layer quasi-geostrophic model by writing the layer 
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interaction matrix L a p(k) as 



L a/3 (k) 



a(k) b(k) 
b{k) a(k) 



(4.4) 



with a(k) and b(k) given by a(k) = k 2 + k 2 R /2 and b(k) — —k 2 R . As will become apparent 
from the argument below, retaining generality and calculating A(fc, q) in terms of a{k) 
and b(k) seems to simplify the details of our calculation. From Eq. ( |3.11[ ) , we write the 
energy spectrum E(k) in terms of the streamfunction spectra Ui(k), Uzft), and C\ 2 {k): 

E(k) = -J2^(k)C a0 {k) (4.5) 

a/3 

= a{k)C u {k) + b(k)C 12 (k) + b{k)C 21 {k) + a{k)C 22 {k) (4.6) 

= a(fc)[t/i(fc) + U 2 (k)} + 2b(k)C l2 (k) (4.7) 

= [a(k) + b(k)]U(k) + b(k)[2C 12 (k) - U(k)}. (4.8) 

Note that we reorganize the expression in terms of 2C\ 2 {k) — U{k) because it is known 
that 2Ci 2 {k) — U(k) < for all wavenumbers k. Likewise, we also know that U(k) > 0. 
Working with spectra that are known to be positive definite or negative definite, makes it 



easier to determine the sign of A(fc, q) in the calculations below. Similarly, from Eq. (3.14), 
the layer- by-layer potential enstrophy spectrum G a {k) reads 

G a (k) = Y, L M k ) L ^(k)C^(k) (4.9) 

Pi 

= L al {k)L al (k)C n {k) + 2L al {k)L a2 (k)C 12 (k) + L a2 (k)L a2 (k)C 22 (k) 7 (4.10) 

and it follows that for the top and bottom layers 

d(fc) = a 2 (k)U 1 (k) + 2a(k)b(k)C 12 (k) + b 2 (k)U 2 (k), (4.11) 
G 2 (k) = b 2 (k)U!(k) + 2a(k)b(k)C 12 (k) + a 2 (k)U 2 (k). (4.12) 

The total potential enstrophy spectrum G(k) = G\{k) + G 2 (k) then reads: 

G(k) = [a 2 (k) + b 2 (k)]U(k) + 4a(k)b{k)C 12 (k). (4.13) 

We use the above equations to evaluate the dissipation rate spectra for the cases of 
streamfunction dissipation and potential vorticity streamfunction. 

4.2. Case 1: Streamfunction dissipation 

We begin with considering the case of streamfunction dissipation in which the top-layer 
has a dissipation operator with spectrum Z?i(fc) = D(k) and the bottom layer has dis- 
sipation operator with spectrum D 2 (k) = D(k) + d(k). For d(k) — we have the case 
of symmetric dissipation in which both layers have the same dissipation operators. For 



d(k) > 0, we have the more general case of asymmetric dissipation. From Eq. (3.201, the 
energy dissipation rate spectrum _D^(fc) is given by 

D E (k) = 2j2D a (k)U a (k) (4.14) 

a 

= 2{D(k)U 1 (k) + [D(k) + d(k)]U 2 (k)}. (4.15) 
This equation can be rewritten in the form 

D E {k)=A^\k)U 1 (k) + A%\k)U 2 {k)+Af{k)[2C l2 {k)-U{k)], (4.16) 
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with Ag\k) and A^\k) given by 

A%\k) = 2D(k) and A% ] {k) = 2[D(k) + d(k)} and A^\k) 
It can also be rewritten in the form 

D E (k) = 4" (k)D(k) + B^ (k)d(k), 

with Bg\k) and Bg\k) given by 



4 1) w = 



2U 1 (k) + 2U 2 (k) = 2U{k) and B^\k) = 2U 2 (k). 



(4.17) 
(4.18) 

(4.19) 



Similarly, from Eq. (3.22|, the potential enstrophy dissipation rate spectrum D G (k) is 
given by 



D G (k) = -2Y^L a 0(k)D (k)C aP (k) 

a/3 







= 2D(k) -J2L a /3(k)C af3 (k) +2d(k) 

a/3 

= 2D(k)E(k) + 2d(k)[b(k)C 21 (k) + a(k)U 2 (k)} 

= 2D{k)E{k) + d(k)b(k)[2C 12 (k) - U(k)} + d(k)b(k)U(k) 

+ 2d(k)a(k)U 2 (k) 
= 2D{k)E(k) + d(k)b(k)[2C 12 (k) - U(k)} + d(fc)6(fc)C/ 1 (/c) 

+ d(k)[2a(k) + b(k)]U 2 (k). 



(4.20) 

(4.21) 
(4.22) 
(4.23) 
(4.24) 



Substituting the energy spectrum E(k) from Eq. (4.8), it immediately follows that D G {k) 
reads 



D G {k) = 4 1) (fc)C/ 1 (fc) + A%\k)U 2 {k) 



AP(k)[2C 12 (k)-U(k)}, 



with A G l \k), A%\k), and A%\k) given by 

A%\k) = 2D{k)[a{k) + b(k)] + d{k)b{k), 

A { (?{k) = 2D(k)[a(k) + b(k)] + d(k)[2a(k) + b(k)}, 

A { G \k) = 2D(k)b{k) + d{k)b{k) = [2D(k) + d{k)]b{k). 



(4.25) 

(4.26) 
(4.27) 
(4.28) 



The potential enstrophy dissipation rate spectrum can also be reorganized in terms of 
D(k) and d(k) and be rewritten as 

D G (k) = B^\k)D(k)+B^\k)d(k), (4.29) 
with B G \k) and B G {k) given by 

B^ ] (k) = 2[a(k) + b(k)]U(k) + 2b(k)[2C 12 (k) - U(k)}, (4.30) 

B^\k) = 2[b(k)C 12 (k) + a(k)U 2 (k)}. (4.31) 

From the above equations we may write the terms of A(fc, q) in terms of the stream- 
function spectra U\(k), U 2 {k) 1 and C\ 2 {k) 1 or in terms of the dissipation operator spectra 
D(k) and d(k). In the former case, we find that A(k,q) reads: 

A(k,q) = k 2 D E (q)-D G (q) (4.32) 
= A^k, qptiq) + A 2 (k, q)U 2 (q) + A 3 (k, q)[2C 12 (q) - U(q)}, (4.33) 
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with Ai(k,q), A 2 (k,q), A 3 (k,q) given by 

A 1 {k,q) = k 2 A^\q)-A%\q) (4.34) 

= k 2 [2D(q)] - 2D(q)[a(q) + b{q)\ - d(q)b(q) (4.35) 

= 2D(q)[k 2 - a(q) - b(q)] - d(q)b(q), (4.36) 

A 2 (k,q) = k 2 A (2) (q)-A {2) (q) (4.37) 

= k 2 2[D(q) + d{q)] - 2D(q)[a(q) + b(q)} - d(q)[2a(q) + b(q)} (4.38) 

= 2D(q) [k 2 - a(q) - b(q)} + d(q) [2k 2 - 2a(q) - 6(g)] (4.39) 

= 2[D(q) + d(q)][k 2 - a(q) - 6(g)] + %)6(g), (4.40) 

A 3 (k,q) = k 2 Af(q)~A^(q) = -[2D{q) + d(q)]b(q). (4.41) 

Reorganizing the terms in the above expression in terms of D(k) and d(k), we can also 
rewrite A(k,q) as 

A(fc, q) = B^k, q)D(q) + B 2 (k, q)d(q), (4.42) 
with B\(k,q), B 2 (k,q) given by 

B 1 (k,q) = k 2 B^(q)-B^(q) (4.43) 

= k 2 [2U(q)} - 2[a(q) + b(q)]U(q) - 2b(q)[2C\ 2 (q) - U(q)} (AAA) 

= 2[k 2 - a(q) - b(q)]U{q) - 2b(q)[2C 12 (q) - U(q)}, (4.45) 

B 2 (k,q) = k 2 B^\q)~B^(q) (4.46) 

= k 2 [2U 2 (q)] - 2[b(q)C 12 {q) + a(q)U 2 (q)] (4.47) 

= 2[k 2 - a{q)]U 2 (q) - 2b(q)C 12 (q). (4.48) 

Using the above expressions for A(fc, g), we will now prove the following statements: 
Proposition [T] establishes the flux inequality for the case of symmetric streamfunction 
dissipation. Proposition [2] gives a sufficient condition in terms of the dissipation operator 
spectra D(k) and d(k) for satisfying the flux equality for the case of asymmetric stream- 
function dissipation. Proposition [3] provides with an alternate sufficient condition for the 
asymmetric case formulated in terms of the streamfunction spectra C\ 2 (q) and U 2 (q). 

Proposition 1. Assume streamfunction dissipation with d(k) = and k 2 — a(q) — 
6(g) < and 6(g) < 0. Then A(k, q) < 0. 

Proof. Under the assumption d(k) = (symmetric dissipation), A(fc, g) reads: 

A(fe,g) = A 1 (fc,g)C/ 1 (g)+^ 2 (fc,g)f/ 2 (g)+A 3 (fc I g)[2C 12 (g)-[/(g)] ) (4.49) 
with Ai(k, q), A 2 (k,q), and A 3 (k, q) given by 

A,(k, q) = 2D(q)[k 2 - a(q) - 6(g)], (4.50) 

A 2 (k, q) = 2D(q)[k 2 - a(q) - 6(g)], (4.51) 

A 3 (k,q) = -2D(q)b(q). (4.52) 

Now we note that D(q) > 0, and by hypothesis, k 2 — a(q) — 6(g) < 0, and therefore we 
immediately obtain that Ax(k,q) < and A 2 (k,q) < 0. We also know, by hypothesis, 
that 6(g) < 0, and therefore we also have A 3 (k, q) > 0. Finally, we know that Ui(q) > 
and U 2 (q) > and 2Ci 2 (g) - U(q) < 0. It follows that all three terms in Eq. ( |4.49[ ) are 
negative, and therefore A(k,q) > 0. □ 
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Proposition 2. Assume streamfunction dissipation with d(k) > and k 2 — a(q) — 
6(g) < and 6(g) < and 

d(q) , k 2 - a{g) - b(q) 
2%) " bjq) ■ (453) 

Then it follows that A(fc, q) < 0. 

Proof. We now write A(fc, q) as 

A(k,q) = A 1 (k,q)U 1 (q)+A 2 (k,q)U 2 (q)+A 3 (k,q)[2C 12 (q)-U(q)}, (4.54) 
with Ai(k, q), A 2 {k,q), and ^(fcjg) given by 

A^q) = 2D(q)[k 2 - a(q) - b(q)} - d(q)b(q), (4.55) 

A 2 (k, q) = 2[D(q) + d(q)][k 2 - a(q) - 6(g)] + d(g)6(g), (4.56) 

A 3 (k,q) = -[2D(q)+d(q)]b(q). (4.57) 

From the assumptions fc 2 — a(g) — 6(g) < and 6(g) < 0, given by hypothesis, and the fact 
that D(q) > and d(q) > 0, we see that A 2 (k, q) < and A 3 (k, q) > 0. Since U 2 (q) > 
and 2Ci 2 (<z) — U(q) < 0, it follows that the second and third terms of Eq. (4.54 1 are 



negative, and therefore A(fe,g) < A\(k, q)U\{q). We also note that U±(q) > and the 



hypothesis given by Eq. (4.53) implies that Ai(k,q) < 0. It follows that A(fe,g) < 0. □ 



Proposition 3. Assume streamfunction dissipation with d(k) > and k 2 — a(q) — 
6(g) < and 6(g) < and C\ 2 {q) < U 2 (q). Then, it follows that A(fc, q) < 0. 

Proof. We write A(k,q) as 

A(fc, g) = B x {k, q)D(q) + B 2 {k, q)d(q), (4.58) 
with Bi(k,q) and B 2 {k,q) given by 

Bt(k, q) - 2[fc 2 - a(q) - b(q)]U(q) - 2b(q)[2C 12 (q) - U(q)}, (4.59) 
B 2 (k, q) = 2[k 2 - a(q)]U 2 (q) - 2b(q)C 12 (q). (4.60) 

From the hypotheses k 2 — a(q) — 6(g) < and 6(g) < 0, and also because U(q) > 
and 2Ci2(g) — U(q) < 0, we immediately find that Bi(k,q) < 0. From the hypothesis 
Ci 2 (q) < U 2 (q), we can also show that 

B 2 (k, q) = 2[k 2 - a(q)]U 2 {q) - 26(g)C 12 (g) (4.61) 

< 2[k 2 - a(q)]U 2 (q) - 2b(q)U 2 (q) (4.62) 

= 2[fc 2 -a(g)-6(g)]C/ 2 (g) < 0. (4.63) 

Since d(q) > and D(q) > 0, it follows that A(fc,g) < 0. □ 

With respect to the hypotheses in the above propositions, we note that for the case 
of the two-layer quasi-geostrophic model, with a(q) = q 2 + k 2 R j2 and 6(g) = — fc^/2, the 
assumptions fc 2 — a(q)— 6(g) = fc 2 — g 2 < and 6(g) < hold mathematically for all k < q. 
The assumption given by Eq. (4.53) in proposition [2] and the assumption C\ 2 {q) < U 2 (q) 



in proposition [3] are the physically substantive assumptions. 

4.3. Case 2: Potential vorticity dissipation 

The case of potential vorticity dissipation is mainly an interesting mathematical curios- 
ity. We are assuming that every layer is being dissipated by a dissipation operator with 
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positive definite spectrum D a (k), applied on the potential vorticity field q a of the corre- 
sponding layer, instead of the streamfunction tp a . The idea here is to determine whether 
this simple change in the configuration of the dissipation terms has any notable impact 
on the robustness of the flux inequality. As we have explained previously, in this case, 
the dissipation matrix D a p{k) is given by D a p{k) = —D a (k)L a p(k), and the energy 
dissipation rate spectrum _Ds(fc) and the potential enstrophy dissipation rate spectrum 



D G (k) are given by Eq. (3.27) and Eq. (3.29) respectively. 

Now, let us restrict ourselves again to the case of a general two-layer model with 



£>i(fc) = D(k) and D 2 (k) = D(k) + d(k). From Eq. ( |3.27[ ), we calculate the energy 
dissipation rate spectrum DE{k) which reads 



D E (k) = -2j2D a (k)L aP (k)C al3 (k) 

a/3 



= 2D(k) 



-J~]L a p(k)C a p(k) 

a/3 

= 2D{k)E{k) + 2d(k)[b(k)C 12 (k) + a(k)U 2 (k) 



(4.64) 



2d(k)[-L 21 (k)C 21 (k) - L 22 (k)C 22 (k)} (4.65) 



(4.66) 



Similarly, from Eq. (3.29), we calculate the potential enstrophy dissipation rate spectrum 

d(k)]G 2 (k) 



DQ(k), which is given by 

D G {k) ^2^D a {k)G a {k) = 2D(k)G 1 (k) + 2[D(k) 

(X 

= 2D{k)[G 1 {k) + G 2 (k)} + 2d(k)G 2 (k). 
Both spectra may therefore be rewritten according to the form 

D E (k) = D(k)B^ (k) + d(k)B^ (k), 
D G (k) = ^(fc)^ (k) + d{k)B^ ] (k), 

with B%\k), B^\k), B£\k), and B^' ' (k) given by: 

B E l \k) = 2E(k) = 2[a(k)U(k) + 2b(k)C 12 (k)} 

= 2[a(k) + b{k)]U{k) + 2b(k)[2C 12 (k) - U(k)}, 



2[b(k)C 12 (k) + a(k)U 2 (k)} 

b(k)[2C 12 (k) - U(k)} + b{k)U{k) + 2a(k)U 2 (k) 

b{k)Ux{k) + [2a(k) + b(k)]U 2 (k) + b(k)[2C 12 (k) 



U(k)}, 



(4.67) 
(4.68) 

(4.69) 
(4.70) 



(4.71) 
(4.72) 

(4.73) 
(4.74) 
(4.75) 



and 



B [ c\k) 



2[G 1 (k)+G 2 (k)} (4.76) 

= 2[a 2 {k) + b 2 (k)]U{k) + 8a(k)b(k)C 12 (k) (4.77) 
= 2[a 2 (k) + b 2 (k)]U(k) + 4a(k)b{k)[2C 12 (k) - U(k)} + Aa{k)b{k)U {k) (4.78) 

= 2[a{k) + b(k)] 2 U(k) + 4a(k)b(k)[2C 12 (k) - U(k)], (4.79) 

B {2) (k) = 2G 2 (k) = 2b 2 (k)U 1 (k) + 4a(k)b(k)C 12 (k) + 2a 2 (k)U 2 (k) (4.80) 
= 2a(k)b(k)[2C 12 (k) - U(k)} + 2a{k)b(k)U (k) + 2b 2 (k)U 1 (k) 

+ 2a 2 {k)U 2 (k) (4.81) 
= 2[a(k)b(k) + 6 2 (fc)]C/i(fc) + 2[a(k)b{k) + a 2 (k)]U 2 (k) 
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+ 2a(k)b(k)[2C 12 (k) - U(k)} (4.82) 
= 2b(k)[a(k) + 6(fc)]C/i(fc) + 2a(k)[a(k) + b(k)]U 2 (k) 

+ 2a(k)b(k)[2C 12 (k) - U(k)}. (4.83) 

Equivalently, the terms of D E {k) and D G {k) can be reorganized in terms of the stream- 
function spectra U\{q), U 2 (q), and 2C\ 2 (q) — U(q), and rewritten as 

D E (k)=A^(k)U 1 (k) + A ( §\k)U 2 (k)+A^(k)[2C 12 (k)-U(k)}, (4.84) 

D G (k) = 4 1) (fc)C/ 1 (fc) + A {2) (k)U 2 (k) + Af{k)[2C l2 {k) - U(k)\, (4.85) 

with A%\k), 4 2) (fc), 4 3) (fc), A%\k), 4 2) (fc), and 4 3) (fc), given by 

4 ) ( fc ) = 2D(k)[a(k) + b(k)] + d(k)b(k), (4.86) 

A {2 \k) = 2D(k)[a(k) + b(k)] + d(k)[2a(k) + b(k)}, (4.87) 

A%\k) = 2D(k)b(k) + d(k)b(k), (4.88) 



and 



A ( e\k) = 2D(k)[a{k) + &(fc)] 2 + 2d{k)b{k)[a{k) + b(k)], (4.89) 

A%\k) = 2D(k)[a(k) + b(k)} 2 + 2d(k)a(k)[a(k) + b(k)], (4.90) 

A^ik) = 4D(k)a(k)b(k) + 2d(k)a(k)b(k). (4.91) 

From the above equations we now write the terms of A(fc,g), collected in terms of the 
streamfunction spectra, as: 

A(k,q) = k 2 D E (q)-D G (q) (4.92) 

= A,(k, qJUtiq) + A 2 (k, q)U 2 (q) + A 3 (k, q)[2C 12 (q) - U(q)}, (4.93) 

with Ai(k, q), A 2 (k,q), A 3 (k,q) given as 

A 1 (k,q) = k 2 A^(q)-A^(q) (4.94) 
= 2k 2 D(q)[a(q) + b(q)} + k 2 d{q)b(q) - 2D{q)[a(q) + 6(g)] 2 

-2d(q)b(q)[a(q) + b(q)} (4.95) 

= 2D(q) [a(q) + 6(g)] [fc 2 - a(q) - b(q)} + d(q)b(q) [k 2 - 2a(q) - 2b{q)] , (4.96) 

A 2 (k,q) = k 2 A^(q)-A^(q) (4.97) 
= 2k 2 D(q)[a(q) + b(q)} + k 2 d(q)[2a(q) + b(q)] - 2D(q)[a(q) + b(q)] 2 

-2d(q)a(q)[a(q) + b(q)} (4.98) 
= 2D(q)[a{q)+b(q)][k 2 - a(q) - b(q)} 

+ d(q){2k 2 a(q) + k 2 b(q) - 2a(q)[a(q) + 6(g)]} (4.99) 
= 2D(q)[a(q)+b(q)][k 2 - a(q) - b(q)} 

+ d(q){2a(q)[k 2 - a(q) - b(q)} + k 2 b(q)} (4.100) 

= [fc 2 - o(g) - b(q)]{2D(q)[a(q) + 6(g)] + 2a(q)d(q)} + k 2 d(q)b(q), (4.101) 

A 3 (fc,g)-fc 2 4 3) (g)-4 3) (g) (4.102) 

= 2k 2 D(q)b(q) + k 2 d(q)b(q) - 4D(q)a(q)b(q) - 2d(q)a(q)b(q) (4.103) 

= 2D(q)b(q)[k 2 - 2a{q)\ + d(g)6(g)[fc 2 - 2a(g)] (4.104) 



the two-layer quasi- geostrophic model 
b(q)[2D(q) + d(q)][k 2 - 2a(q)}. 
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(4.105) 



Similarly to the previous case of the streamfunction dissipation, we now derive the 
following two propositions. Proposition [4] establishes the flux equality for the case of 
symmetric potential vorticity dissipation. Proposition [5] provides a sufficient condition for 
satisfying the flux inequality, for the case of asymmetric potential vorticity dissipation. 

Proposition 4. Assume potential vorticity dissipation with d{k) — for all k > and 
a(q) + b(q) > and k 2 - a(q) - b(q) < and b(q) < 0. Then, it follows that A(k, q) < 0. 

Proof. Under the assumption that d{k) — for all wavenumbers k, the general form 
of A(k,q) and the coefficients Ai(k,q), A 2 (k,q), A 3 (k,q) simplify to: 

A(k,q)=A 1 (k,q)U 1 (q)+A 2 (k,q)U 2 (q)+A 3 (k,q)[2C 12 (q)-U(q)], (4.106) 

A x {k,q) = 2D(q)[a(q)+b(q)}[k 2 - a(q) - b(q)}, (4.107) 

A 2 (k,q) = 2D(q)[a{q) + b(q)][k 2 - a(q) - 6(g)], (4.108) 

A 3 (k,q) = 2D{q)b{q)[k 2 - 2a{q)\. (4.109) 

From the assumptions a(q) + b(q) > and k 2 — a(q) — b(q) < 0, it immediately follows 
that Ai(k,q) < , and A 2 (k, q) < 0. We also find that a(q) > —b(q) > 0. Using this 
result in conjugation with the previous two assumptions, we show that 



2a(q) = [k 2 - a(q) - b(q)} - a(q) 
< [k 2 - a(q) - b(q)} + b(q) 
<k 2 - a(q) - b(q) < 0. 



b(q) 



(4.110) 
(4.111) 
(4.112) 



It follows that b(q)(k 2 - 2a(q)) > 0, and therefore A 3 (k,q) > 0. Since Ui(q) > and 



U 2 {q) > and 2C\ 2 {q) — U{q) < 0, it follows that all 3 terms in Eq. (4.1061 are negative 
and therefore A(k,q) < 0. □ 

Proposition 5. Assume potential vorticity dissipation and assume that a(q)+b(q) > 
and k 2 — a(q) — b(q) < and b(q) < 0. We also assume that d(q) satisfies 

[a(q) + b(q)}[k 2 -a(q)-b(q)} 



d(q) 



< 



D(q) 

Then, it follows that A(k,q) < 0. 



b(q)[k 2 - 2a(q)-2b(q)} 



(4.113) 



Proof. Recall that the general form of A(k,q) reads 

A(k,q)=A 1 (k,q)U 1 (q) + A 2 (k,q)U 2 (q)+A 3 (k,q)[2C 12 (q)-U(q)}, (4.114) 

with Ax(k, q), A 2 (k,q), A 3 (k,q) given by 

A 1 (k,q) - 2D(q)[a(q) + b{q)][k 2 - a(q) - 6(g)] + d{q)b(q)[k 2 - 2a(q) - 2b(q)}, (4.115) 

A 2 (k, q) = [k 2 - a(q) - b(q)]{2D(q)[a(q) + 6(g)] + 2a(q)d(q)} + k 2 d{q)b{q), (4.116) 

A 3 (k,q) = b(q)[2D(q) + d{q)][k 2 - 2a(q)}. (4.117) 

As in the previous proof, our assumptions also imply that a(q) > and k 2 — 2a(q) < 0. 
These are sufficient for determining the sign of A 2 (k, q) and A 3 (k, q). For the case of 
A 2 (k,q), we note that it consists of two terms, both of which positive, as the product 
of a negative and a positive factor. It follows that A 2 (k,q) < 0. Likewise, A 3 (k,q) is 
the product of two negative and one positive factors, and therefore A 3 (k, q) > 0. Since 
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U 2 {q) > and 2C 12 (q) - U(q) < 0, it follows that A(k,q) < A^qjU^q). For the 
case Ai(k,q), we see that it consists of two competing terms: the first negative and the 
second positive. However, from Eq. (4.113), it follows that A\[k,q) < 0, consequently 
A(k,q) < 0, as claimed. □ 

Comparing the above propositions with propositions [l] [2j and |3j corresponding to the 
case of streamfunction dissipation, we see that the additional assumption a(q) + b(q) > 
is needed. For the case of the two-layer quasi-geostrophic model, we have a(q) + b(q) = 
q 2 > 0, so the assumption is satisfied. We also note that Eq. (4.53) and Eq. (4.113) 
are not identical, but they are nonetheless very similar. Finally, I was not able to find 
an appropriate counterpart to proposition [3| Whether there exists such a counterpart 
remains an open question. 



5. Summary of main results 

In the previous section we derived a series of propositions that provide us with sufficient 
conditions for establishing the inequality A(fc, q) < for all wavenumbers q > k. These 
results are mathematically rigorous and require no physical assumptions. As explained in 
the beginning, this inequality in turn establishes the flux inequality k 2 HE(k) — Hdk) < 
0. In this last step we have to assume that the wavenumber k is not in the forcing 
range. In quasi-geostrophic turbulence models, the forcing spectrum is controlled by the 
baroclinicity instability. Consequently, implicit in the argument below is the assumption 
that baroclinic instability is negligible at large wavenumbers. This has been originally 
proposed by Salmon (1978 |1980 ) and has been corroborated by Welch & Tung (1998) 
and Tung & Orlando (2003). For investigative purposes, in a numerical simulation, we 
can assume control of the forcing spectrum by using antisymmetric random gaussian 



forcing, as explained by Gkioulekas (20121 



Our results are restricted to the two-layer model under two distinct dissipation config- 
urations: streamfunction dissipation and potential vorticity dissipation. For both config- 
urations, we have shown that as long as the same dissipation operator is applied on both 
layers, the flux inequality is satisfied for all wavenumbers not in the forcing range. This 
follows from propositions [T] and |4j For the case of asymmetric dissipation, in which the 
operators for the two layers are not identical, we need to distinguish between the case of 
streamfunction dissipation and potential vorticity dissipation. 

Under asymmetric streamfunction dissipation, in which the dissipation operator of the 
top layer has spectrum D\(k) — D(k) and the dissipation operator of the bottom layer 
has spectrum D 2 (k) = D(k) + d(k), we have shown by proposition [2] that 



d(q) < k 2 - a(q) - b(q) 



2D(q) 



Kq) 



A(k,q)<0. 



(5.1) 



For the case of the two-layer model, using a(q) = q 2 + k 2 R /2 and b(q) = —k R /2, gives 
k 2 — a(q) — b(q) = k 2 — q 2 and consequently the above statement reduces to 



d(q) 4(g 2 - k 2 ) 
D(q) 



A(k,q)<0. 



(5.2) 



For the typical case in which both layers are dissipated by the same hyperdiffusion 
operator, and the bottom layer is dissipated by an Ekman term, we have D(q) = vq 
and d(q) — VEq 2 which gives the following statement 



2p+2 



k\v E < Avq 2 P{q 



A(k,q)<0. 



(5.3) 
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For finite Ekman coefficient ve, it is easy to see that the inequality is violated when k and 
q are too close to each other, so in principle, A(fc, q) could become positive for k ~ q and 
then transition to negative for k <C q- On the other hand, if the wavenumber k is deep 



within the inertial range, the contributions of q ~ k to the integral in Eq. (4.3 1 are bound 
to be negligible since dissipation is not expected to be dominant in the inertial range. 
Furthermore, the inequality only ensures the negativity of an upper bound of A(k,q), 
namely the contribution of the Ai(k, q)Ui(q) term, which is the only term that can be 
positive or negative. Even if the A\(k, q)U\(q) term is positive, for small enough Ekman 
coefficient Ve, it is reasonable to expect it to be overtaken by the other two negative 
terms, since for ve = 0, the Ai(k, q)U\(q) term is also negative. Consequently, we believe 
that for small enough ve, the flux inequality will continue to hold. 

This is as much as can be said without making any phenomenological assumptions. If 
we had a more detailed knowledge of the phenomenology of the streamfunction spectra 
Ui(q), U2(q), and Ci2(q), it would be possible to calculate A(k, q) explicitly and perform 
a more precise investigation. In connection with this matter of the phenomenology of 
the streamfunction spectra, it is worth noting that we have also shown via proposition [3] 
that another sufficient condition for establishing the flux inequality at wavenumber k is 
to have C\2{q) < £^2(9) f° r an wavenumbers q > k. 

For the case of potential vorticity dissipation, from proposition [5] we have a similar 
result that 

d(g) -[a(q) + b(q)} [k 2 - a(q) - b(q)} 

W) - Kq W -2a { q)-2 m => A(fc ' 9) * °< (5 ' 4) 

and setting a(q) and b(q) to a(q) = q 2 + k R /2 and b(q) — — fcJj/2, as required for the 
case of the two-layer model, gives a(q) + b(q) = q 2 , k 2 — a(q) — b(q) = k 2 — q 2 , and 
k 2 — 2a(q) — 2b(q) = k 2 — 2q 2 , and the general result reduces to 

d(q) ^ A(q 2 -_k 2 ) q 2 
<7? 



m- ^ 2q 2 -k 2 =^ A ^^ - ^ 



A typical choice for the dissipation operator spectra D(k) and d(k) is to choose D(q) = 
vq 2p and d(q) = ve- This choice gives the statement 

v E k 2 R (2q 2 - k 2 ) < Avq 2p+2 {q 2 - k 2 ) => A(k, q) < 0. (5.6) 

Similar considerations apply to this case, as in our discussion above of the case of stream- 
function dissipation. However, the sufficient condition here gives a tighter inequality, 
which suggests that a violation of the flux inequality may be easier under potential 
vorticity dissipation. 

It is interesting to note that if one uses differential hyperdiffusion, whereby the hyper- 
dissipation term at the bottom layer has greater hyperviscosity coefficient than the hy- 
perdissipation term at the top layer, then for D(k) = vq 2p+2 and d(k) = VEq 2 + Avq 2p+2 , 
we have, for the case of streamfunction dissipation, the statement 

k 2 R v E < q 2p [4v(q 2 - k 2 ) - k 2 R Av] A(fc, q) < 0, (5.7) 

and, for the case of potential vorticity dissipation, using D(k) = vq 2p and d(k) — ve + 
Avq 2p , the statement 

v E k 2 R (2q 2 - k 2 ) < q 2p [{q 2 - k 2 ){Avq 2 - Avk 2 R ) - 2k 2 R Avq 2 \ A(fc, q) < 0. (5.8) 

Here, we have assumed that the hyperviscosity coefficient at the top layer is v and at the 
bottom layer is v + Av with Av > 0. We see that in both cases, the inequality constraint 
on ve becomes tighter. We may therefore speculate that with increasing Av, it may 



20 



Eleftherios Gkioulekas 



become easier to break the flux inequality in the two-layer quasi-geostrophic model. This 
possibility warrants further numerical investigation. 



6. Conclusion and Discussion 

The flux inequality and the possibility of its violation was discussed previously in a 
less systematic manner by Gkioulekas & Tung (2007a). In the present paper, we have 



provided a more systematic treatment of the energy and potential cnstrophy dissipation 
rate spectra for the general case of an arbitrary multi-layer model. We gave a more 
careful derivation of the sufficient conditions for satisfying the flux inequality, for the 
case of streamfunction dissipation, and extended these results to the case of potential 
vorticity dissipation. The careful form of our argument makes it possible to consider 
the possibility of differential hyperdiffusion, and we now have an ideal launching point 
for a more thorough investigation of the flux inequality. Our results on the non-trivial 
dependence of the dissipation rate spectra of energy and potential enstrophy on the 
energy and potential enstrophy spectra via the streamfunction spectra, are also very 
relevant to the correct formulation of closure models for multi-layer quasi-geostrophic 
systems. We would like to stress again that these results are mathematically rigorous and 
the only underlying assumption is that forcing, which is driven by baroclinic instability, 
is confined to large scales. This assumption is not needed to establish propositions 1-5, 
but it needs to be introduced when taking the next step to establish the flux inequality 
from the conclusion of these propositions. No phenomenological assumptions about any 
spectrum are needed at any step of the argument. 

One limitation of the current investigation is that we have disregarded the beta term, 
mainly to avoid the mathematical difficulties associated with the anisotropic nature of 
the term. This elimination can be tolerated, from a physical standpoint, as long as the 
beta term is active only in the forcing range and the baroclinic forcing at the same forcing 
range is powerful enough to overshadow the beta term. If barotropization of the two-layer 
quasi-geostrophic model is indeed catalyzed by the beta term, as suggested by |Venaille,1 
Vallis & Griffies (2012), then the effect of the beta term will probably be to strengthen 
the robustness of the flux inequality, with increasing values of the /3 coefficient. This 
tendency can be resisted by asymmetric dissipation, but it is not obvious whether the 
two tendencies can counterbalance each other for any value of the (3 coefficient, or which 
of the two tendencies gets to dominate for realistic choices of (3. On the other hand, as 
long as the effect of the beta term remains limited to large scales (i.e. planetary and 



synoptic scales), it will not contribute to the integrals of Eq. (4.1) and Eq. (4.2) and the 
results reported in this paper will remain entirely unaffected. 

As was explained previously in the introduction of the present paper, the |Tun g & 
Orlandol (120031) simulation indicates that the downscale energy cascade can indeed be 



uncovered. However, in the absence of the theory developed in the present paper and 
previous papers |Gkioulekas||2012| |Gkioulekas fc Tung||2005a|6[ |2006| |2007a|6[ ), |Tung 



& Orlando (2003) were not motivated to explore the layer-by-layer phenomenology. A 



numerical investigation of this phenomenology would go a long way towards building 
our understanding of the dynamics of the two-layer quasi-geostrophic model, as a step- 
ping stone for further investigation of multi-layer quasi-geostrophic systems. Just as 
two-dimensional turbulence research turned out to be interesting, exciting, and full of 
surprises, in its own right, we believe that quasi-geostrophic models, such as the two-layer 
model, will prove to have an even richer set of surprises in store for us. The two-layer 
model is a very good platform for investigating the dynamics of coexisting cascades, and 
for that reason alone, further investigation is warranted. 
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The idea of a flux inequality was first brought up in email communication between 
Sergey Danilov with the author and Ka-Kit Tung, in the context of two-dimensional 
Navier-Stokes turbulence. 



Appendix A. Proof that E(k) is positive definite 

In this appendix we prove the inequality 2\C a p(k)\ < U a {k) + Up(k), and we use it to 
show that if the matrix L a p{k) satisfies the following conditions: 

L a p(k) > 0, for (Al) 

2L aa (k) + J2(LMk)+L Pa (k))<0, (A 2) 

a/3 

then the energy spectrum E(k) will be positive definite with E(k) > 0. 
To establish the inequality, we first note that 

(ip a ±ipp,'ip a ±il)p) k = (ip a ,ip a ) k ± 2 (ip a ,ipp) k + (il>p,il>p) k (A3) 
= U a (k) + Up{k)±2C a p(k). (A 4) 

Since (tp a ±ipp,ip a ± ipp) k is by definition positive definite, it follows that U a (k)+Up(k) — 
2C a p(k) > and U a (k) + Up(k) + 2C a p(k) > 0, and therefore we have 

2\C a p{k)\ <U a (k) + Up(k). (A 5) 

Now, to show that the energy spectrum E(k) is positive definite, we begin by rewriting 
Eq. pill as follows: 



2E(k) = -2j2L a p(k)C a p(k) = -2j2L aa (k)U a (k) - 2 £ L a p{k)C a p{k) (A 6) 

a/3 a a/3 

a^/3 

= -2j2L aa (k)U a {k) - L a p{k)p a {k) + Up(k)} (A 7) 

a a/3 
ajtp 

- J2 L a p (k) [2C a p (fc) - U a (fc) - Up (k)\ (A 8) 

a/3 
a#/3 

= -2j2L aa (k)U a {k)-J2[L a p(k)+Lp a (k)]U a {k) (A9) 

a a/3 

- J2 L a p{k)[2C a p(k) - U a {k) - Up(k)\ (A 10) 

a/3 
a//3 

= -^{2L aQ (fc) + ^\L a p{k) + Lp a (k)]}U a (k) (All) 

a 8 
a#/3 

- J2 L a p(k)[2C a p{k) - U a (k) - Up(k)\. (A 12) 

a/3 
a#/3 

We note that U a (k) > 0, since U a (k) is positive definite, and we have just shown that 
2C a p(k) — U a (k) — Up(k) < 0. Combining these with the assumptions given by Eq. (All 
and Eq. (A2), we see that both terms in our expression for 2E(k) are positive and 
therefore E(k) > 0. 
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For the case of a two-layer model with general matrix L a p(k) given by 



(A 13) 



the conditions given by Eq. (Al) and Eq. (A2) reduce to b(k) < and a(k) + b(k) > 0. 
For a two-layer quasi-geostrophic model, we have a(k) = k 2 + k 2 R /2 and b(k) = —k 2 R /2 
and both conditions are readily satisfied. 



Appendix B. Derivation of dissipation rate spectra 

In this appendix, we will show that the energy dissipation rate spectrum DE{k) and 
the layer-by-layer potential enstrophy dissipation rate spectra Dq o (k) are given by 



D E (k) =2j2D a p(k)C a0 (k), 

uB 



{k)D a ^[k)C M {k). 

Hi 



(Bl) 
(B2) 



The proof mirrors the argument used by Gkioulekas (2012) to derive the energy forcing 



spectrum and the potential enstrophy forcing spectrum for the same model. We begin 
by writing the governing equation for the streamfunction field ipa as 

9 4r + E KpWei 9p) = E Kftihlh + E *$h- 

8 MP 
Differentiating the streamfunction spectrum C a p(k) with respect to time gives 



(B3) 



and we may write a governing equation for C a g(k) in the form: 



dt 



(B4) 



(B5) 



Here, T Qj g(fc) is the contribution from the nonlinear Jacobian term, D a p(k) is the con- 
tribution from the dissipation term, and $ a p{k) is the contribution from the forcing 
term. The dissipation term T) a p{k) can now be obtained by replacing in Eq. (B4) the 
streamfunction time-derivative dtpa/dt with the dissipation term -y-^aB ^M^v This 
gives 

^ (*) = - ( E ^1 %S^S , i>? ) - ( Tpa , E ) (B 6) 

= -Z)[^(*)^( fc ) C ^( fc ) + ^7 1 ( fe ) £, 7«(*) Cf a«(fc)]. (B7) 
7(5 

We may now easily write the dissipation rate spectra and D(j(k) by applying on 

D a p(k) the linear operators indicated by Eq. (3.111 and Eq. (3.14). We therefore find 



that the energy dissipation rate energy spectrum Z?e(fc) is given by 
D E {k) = - WO^M*) 

a/3 



(B» 
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= ^[L Q/3 (fc)L- 7 1 (fc)^(fc)C^(A : )+i^(A : )L^ 1 (A : )^ 75 (A : )C ct5 (fc)] (B9) 

= S Pl D lS (k)C ps (k) + 5 ai D 1& {k)C a5 {k) (B 10) 

= D ps {k)C ps (k) + D l5 (k)C lS {k) = 2j2D a0 (k)C aP (k). (B 11) 

06 7(5 a/3 

The layer-by-layer potential enstrophy spectrum Dc a (k) is likewise given by 

D Ga (k) = L aP (k)L ai (k)'D M (k) (B 12) 

/3 7 

= -J2 L a p(k)L aj (k)[L^(k)D Se (k)C 7S (k) + L-l(k)D Se {k)Cp £ {k)] (B13) 

fJfSe 

= -J2 5 aS L ai (k)D S£ {k)C ie (k) - ]T 5 a sL a p(k)D Se (k)Cp e (k) (B 14) 

•ySe fi5e 

= -J2 L ai {k)D ae (k)C 1£ {k) - ]T L afs (k)D ae (k)C fjE (k) (B 15) 

7e fie 

= -2^L a0 (k)D ai (k)C 0J (k). (B16) 

/3 7 

The corresponding conservation laws read 

«f> + ««*> --!>,(*) +»(»). (BIT) 

Mfi + ^ = _ DcW + fc(t) . (B18 ) 

We see that positive Dsik) and Da(k) correspond to the case where the dissipation 
terms are truly dissipative. This concludes the argument. 
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